Challenges in continuum modeling of intergranular fracture 
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Intergranular fracture in polycrystals is of- 
ten simulated by finite elements coupled to a 
cohesive-zone model for the interfaces, requiring 
cohesive laws for grain boundaries as a function 
of their geometry. We discuss three challenges 
in understanding intergranular fracture in poly- 
crystals. First, 3D grain boundary geometries 
comprise a five dimensional space. Second, the 
energy and peak stress of grain boundaries have 
singularities for all commensurate grain bound- 
aries, especially those with short repeat distances. 
Thirdly, fracture nucleation and growth depends 
not only upon the properties of grain boundaries, 
but in crucial ways on edges, corners, and triple 
junctions of even greater geometrical complexity. 
To address the first two challenges, we explore 
the physical underpinnings for creating functional 
forms to capture the heirarchical commensurabil- 
ity structure in the grain boundary properties. 
To address the last challenge, we demonstrate a 
method for atomistically extracting the fracture 
properties of geometrically complex local regions 
on the fly from within a finite element simulation. 

The nucleation and propagation of cracks in practical 
engineering materials depends strongly on the mesoscopic 
structure; grain boundaries, polyphase inclusions, dislo- 
cations and other defects determine the toughness. Can 
continuum computational modeling be used to quantita- 
tively study such complex failure modes? 

Consider brittle intergranular fracture — rupture at the 
boundary between two crystallites. Ignore for the mo- 
ment problems like embrittlement caused by impurity 
segregation to grain boundaries, and assume a clean, 
single-phase, equilibrium grain boundary. Direct atom- 
istic simulations are infeasible for anything larger than 
nanocrystals; even simulations that focus on the bound- 
aries [1] will be overwhelmed by the number of relevant 
atoms for systems larger than microns in scale. Hence 



let us imagine a finite-element simulation of the poly- 
crystal coupled to, say, cohesive-zone models for each 
interface 0, H m- (A cohesive law gives the crack 
opening as a function of the traction across the interface; 
it is often parameterized by a peak stress and a total 
energy associated with cleaving [7[.) Can one use atom- 
istic models to measure the fracture properties of the 
individual boundaries, and then use these properties in a 
realistic continuum fracture simulation? We outline here 
some serious challenges involved in continuum modeling 
of intergranular fracture; complete details are available 
in longer publications @, Q . 

The first challenge is that of geometrical complexity. 
The cohesive law will depend on the structure of the 
grain boundary. The macroscopic geometry of a 3D grain 
boundary depends on five parameters that describe the 
relative orientations of the two grains. (The fracture dy- 
namics may in principle depend on properties of the crack 
that are not treated explicitly by the cohesive zone model, 
such as the orientation of the crack front within the grain 
boundary or the three separate stress intensity factors.) 
The atomistic structure also depends on how the two 
crystal lattices are translated with respect to one another 
along the three directions, which can greatly affect the 
pattern of atoms along the boundary and hence the peak 
stress and energy Q . One particular shift will constitute 
a global energy minimum corresponding to the most nat- 
ural configuration Q. 

In a polycrystal, one grain will have to find an energy 
minimizing configuration with several other, neighboring 
grains. For a particular grain boundary in a polycrystal, 
where each grain has been pinned by other neighboring 
grains, there will be a competition between elastic strain- 
ing and plastic deforming. For thick grains, in equilib- 
rium, and away from intersections, one can show that it 
is advantageous for the crystallites to strain slightly to 
allow the boundary to find the global energy minimum. 

To compute the cohesive properties of grain bound- 
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FIG. 1: Commensurate 2D Grain Boundaries. The set 

of points above represent all 2D grain boundary geometries 
that can be simulated in a periodic box of 70 lattice constants 
or less, with a strain of 0.05% or less. 6\ and O2 are the 
tilt angles that define the grain boundary geometry. There 
are gaps near perfect crystals, symmetric grain boundaries, 
and high symmetry grain boundaries because creating a new, 
nearby geometry requires adding flaws (extra dislocations) at 
large separations. [6[ 

aries efficiently, it is useful to use periodic boundary con- 
ditions in directions perpendicular to the grain bound- 
ary, which demands that the the two crystals have fi- 
nite repeat distances along the interface, and that the 
repeat distances be commensurate with one another. We 
have found a systematic method of finding commensurate 
grain boundaries |6[ , and have also generalized it to allow 
for slight elastic strains to mesh the two crystal bound- 
aries together. We can approximate commensurate grain 
boundaries by allowing small strains in either direction. 
Fig. [T] describes the commensurate grain boundaries for 
2D hexagonal crystals; Fig. [2] shows a cross section of the 
five dimensional space of commensurate and near com- 
mensurate grain boundaries for three-dimensional FCC 
crystals. 

These commensurability questions are not only of prac- 
tical importance in efficiently computing the properties 
of grain boundaries; commensurate grain boundaries (es- 
pecially those with short repeat distances) also have espe- 
cially low energies and high peak stresses 0, Q . Modeling 
the geometry dependence of the peak stress and energy 
could be relatively straightforward if they depended in 
comprehensible ways on the five geometrical parameters. 
The second challenge in continuum models of fracture is 
to incorporate the singularities associated with commen- 
surate geometries into appropriate functional forms. 

It is well established that cusp singularities in the en- 
ergy occur at special high- symmetry grain boundaries 
with low repeat distances 0, @, 0, [HI- The cusps 
in energy can be understood by thinking of a high- 
symmetry boundary as an undeformed reference crys- 
tal @. Nearby grain boundaries (whose crystallites 
are rotated by a small angle from the high-symmetry 
boundary) are thus described by decorating the high- 
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FIG. 2: Commensurate 3D Grain Boundaries. The 

set of points above represent a stereographic projection (de- 
fined for normal vector ft by ft = (2, j, k)/ \J i 2 + j 2 + k 2 = 
(sin a cos /3, sin a sin /3, cos a)) of all FCC surface orientations 
that are commensurate with the [100] surface within an area 
less than 100 square lattice constants with a strain of .1% or 
less. This constitutes a three dimensional cross section of the 
five dimensional space of grain boundary geometries, where 
the twist angle parameter has been collapsed onto a 2D plot. 

symmetry boundary with a few extra dislocations, just as 
a low-angle grain boundary in a crystal can be described 
as an array of well-separated dislocations. This analogy 
leads to a functional form for grain boundary energy as 
a function of tilt angle in which the cusps around the 
special high symmetry grain boundaries have the same 
OlogO form as low angle grain boundaries [6]. Fig. |3(a)| 
shows the results for a systematic study of symmetric 2D 
grain boundaries and the resulting fitting function @. 

For the peak stress, it is known that there are jumps at 
the same special grain boundaries [6, 8, 10]. We can also 
understand these jumps by using the dislocation picture 
described above [6[ . As we add a dislocation to the high 
symmetry grain boundary, we add a nucleation point for 
fracture, causing a discontinuity in the peak stress. As 
a result, the plot of peak stress vs. tilt angles is dis- 
continuous at every commensurate geometry (Fig. |3(b)| ). 
By considering the elastic interaction between the extra 
dislocations, we have been able to understand also the 
dependence of the peak stress in the vicinity of the high- 
symmetry boundaries (see Fig. |3(b)] and p). 

Are cohesive laws enough? We have studied this ques- 
tion computationally [7] by comparing a direct atom- 
istic simulation of polycrystalline fracture with a finite- 
element simulation of the same geometry using cohesive- 
law parameters derived from the same interatomic po- 
tential. Fig. [4] shows a snapshot of the two simulations of 
polycrystalline fracture in Stillinger- Weber silicon. Both 
in this case and for other simulations, the atomistic sim- 
ulations fail at significantly lower stresses than the con- 
tinuum simulations. Crack nucleation in both atomistic 
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FIG. 3: Singularities at Special Grain Boundaries. 

Cusps in energy appear at high symmetry grain boundaries 
(figure a) and have the same log shape as the energy of low 
angle grain boundaries. The red line is the functional form de- 
scribed in [6(| . The peak stress as a function of tilt angle (figure 
b) is discontinuous everywhere, with higher values at special 
tilt angles corresponding to high symmetry grain boundary 
geometries. The dependence of peak stress on angle near the 
high symmetry grain boundaries (lines and parabolas shown) 
can be explained using the interactions of the extra disloca- 
tions added |6] . The inset shows the peak stress for the grain 
boundary with tilt angle 49.1 and the nearby geometries. 



and continuum simulations happens not in the middle 
of grain boundaries, but at triple junction lines, edges, 
and corners which are not quantitatively described by 
the cohesive laws for the grain boundary interfaces. Sim- 
ilarly, quantitative understanding of how the crack turns, 
branches, or goes intragranular (Fig. 2]) at triple junctions 
demands that we understand the effects of the irregular 
atomistic configurations at these junctions. 

The third challenge is thus to develop an effective com- 
putational method for modeling more complex local ge- 
ometries. One in principle could incorporate an analyti- 
cal understanding of these local geometries into, for ex- 
ample, cohesive laws for triple junctions, but the geomet- 
rical complexity would seem daunting. Can one rely here 
on direct atomistic simulations? Brittle crack nucleation 




(a) Atomistic simulation at 
10% strain 



(b)FEM simulation at 13.1% 
strain 
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FIG. 4: Comparison of Atomistic and Continuum 
Fracture Simulations. The figures above show a center 
cross-section of the atomistic (figure a) and continuum (figure 
b) poly crystal simulations. The geometry is that of a cube in 
cube with three crystal orientations for the upper half of the 
outer cube, lower half of the outer cube, and the inner cube. 
An upwards displacement is imposed on the top face. The 
color scale indicates the vertical component of stress. In the 
FEM simulation (figure b), we include an interface through 
the center of the inner cube to allow for intragranular frac- 
ture. The atomistic simulation fails at a strain 0.02 smaller 
than the FEM simulation (a percentage error of 20%, with the 
crack propagating partially intragranularly through the inner 
cube and partially through the interface at the top of the in- 
ner cube. The continuum simulation cracks straight through 
the center plane of the inner cube. [7[ 



is a local phenomenon, and the intersection of a growing 
crack with a triple junction edge again will generically 
happen at a point. A feasible atomistic simulation of 
the local region of interest could be launched whenever 
the continuum simulation reached a stress state where 
its cohesive laws become unreliable. The information 
about the local geometry (elastic strains, grain orienta- 
tions, and impinging crack surfaces) would be transferred 
from the continuum simulation to generate the atomistic 
configuration, and the results of the atomistic simula- 
tion (nucleation thresholds, crack branching and turning 
events) passed back to the finite-element simulation. 

As an example of such a method, Fig. [5] shows an 
Overlapping Finite Elements and Molecular Dynamics 
(OFEMD) simulation of fracture at a triple-junction, 
generated automatically in this fashion p}. The grains 
were generated using geometry, orientation, and bound- 
ary strains passed from a finite element simulation, run- 
ning on a separate machine and communicating either 
through a command-line interface, a Web service, or a 
database. OFEMD uses the DigitalMaterial [l2[ atom- 
istic simulation environment to deform and relax the 




(a) Initial configuration (b) Cracked configuration 



FIG. 5: Direct Atomistic Simulation of Triple Junc- 
tions. A cylinder of atoms decorates a finite element edge 
that coincides with a triple junction of grains. Figure a shows 
the initial configuration. Figure b shows the atoms after they 
have been deformed according to the finite element displace- 
ments and relaxed, opening a crack along two of the interfaces. 

atomistic coordinates, allowing the failure information to 
be recorded. OFEMD can be downloaded from p~3( ] . 

We have discussed three main challenges involved in 
continuum modeling of polycrystal fracture. First, ex- 
ploring the cohesive properties of 3D grain boundaries 
involves exploring a 5D space. Second, the peak stress 
and energy have singularities at all commensurate grain 
boundaries. Even if it weren't for the first two challenges, 
our comparisons of atomistic and finite element simula- 
tions of polycrystal fracture show that cohesive proper- 
ties of the interfaces alone are not enough to model the 
fracture of polycrystals using continuum methods. Sites 
such as triple junctions, edges, and corners of grains are 
important nucleation sites. In order to resolve this last 
challenge, we suggest the use of direct atomistic modeling 
of local regions of interest. 
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